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Abstract 

We apply Kauffman's automata on small-world networks to study the crossover be- 
tween the short-range and the infinite-range case. We perform accurate calculations 
on square lattices to obtain both critical exponents and fractal dimensions. Partic- 
ularly, we find an increase of the damage propagation and a decrease in the fractal 
dimensions when adding long-range connections. 
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1 Introduction 



The Kauffman model [1] or more generally random Boolean networks have 
been studied in the past to describe genetic regulatory networks but are in fact 
very general because they do not assume any particular function of the nodes 
and can be applied on any topology. While in the late eighties much work has 
been done on the damage spreading transition renewed interest in this model 
has risen in the last years because of the issues of sychronisation[2], stabil- 
ity [3], control of chaos [4] and investigations on scale-free topologies [5, 6]. A 
topology that is very characteristic for human networks is the small-world as 
introduced by Watts and Strogatz [7] . It is the purpose of this paper to present 
some numerical studies of the phase transition on the square lattice, comple- 
menting work by Stauffer [8], and then on its small- world variant showing the 
appearance of a new universality class. 
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Let us first review Kauffman cellular automata. The Kauffman model is a 
random mixture of all possible Booleaan rules. Each of N lattice sites hosts 
a Boolean variable ( spin up or down) which is either zero or unity. The 
time evolution of this model is determined by N functions /j (rules) which are 
randomly chosen for each site independently, and by the choice of K input 
sites {jxii)} for each site i. Thus the value <7j at site i for time t+1 is given 
by: 

a, t (t + 1) = fi{a h {t), a jK {t)) (i = 1, 2, . . . , N) (1) 
Each Boolean function /j is specified, once its value is given for each of the 
2 K possible neighbour configurations. Here we are concerned with those cases 
where both the inputs and the chosen Boolean functions do not change with 
time (quenched model). 

Regular Small-World Random 




The Kauffman model has usually been considered in the last years, on one 
hand with an infinite range of interaction, known as the 'mean-field' approx- 
imation [9,10] in which for each site i , the K inputs are randomly chosen 
among the N spins. Such approximation was introduced to describe genetic 
interactions, to explain mutations and stability processes in biological genes. 
Indeed, this was the original proposal of the Kauffman model. On the other 
hand, this model has also been simulated on lattices with nearest-neighbour 
interactions, the so-called short-range case [8]. If one considers a square lattice, 
as we do here, for each site % the K inputs are just the four nearest neighbours, 
thus there being 2 4 = 16 neighbor configurations, and therefore 2 16 = 65536 
different possible rules. One can implement this by introducing a probability p 
of selecting a rule with the result spin up for each neighbor configuration and 
each lattice site. Of course, because of the symmetry of Boolean functions, 
one also has a probability 1-p that the result is spin down. 

In the past, the connection topology had been assumed to be either com- 
pletely random or completely regular. But many biological, technological and 
social networks lie somewhere between these two extremes. The 'small-world' 
topology seems to be very interesting for this purpose. In order to get more 
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knowledge about this intermediate range, we interpolate between regular and 
random networks by starting from a square lattice with N sites and K=4 
inputs of nearest-neighbors and introducing a rewiring probability q to each 
input K connecting one site to any other, as indicated in Fig.l. Dynamical 
systems with small-world coupling exhibit enhanced signal-propagation speed 
and synchronizability. In particular, small damages, for example mutations on 
genetic material, spread more easily on small-world networks than on regular 
lattice. 

In the present work, we have performed simulations based on a parallel bit 
manipulation technique called multi-spin coding [11] to attempt to determine 
the critical point p c in analogy to percolation, both for the short-range case 
(q = 0) and for the small-world case (q ^ 0). To this end, we plot the nor- 
malized damage mass (see next section), i.e. our order paramenter ip , versus 
the parameter p. Once determinated the critical points p c for each case, we 
evaluate the critical exponents and the fractal dimensions. 



2 'Damage Spreading' on the Square Lattice 




Fig. 2. Plot of ij) vs p for the short-range 
case (q = 0.00). Lattice sizes L=256; o, 
L=384; A and L=512; O. 




Fig. 3. Plot of ip vs p for the small-world 
case (q = 0.05). Lattice sizes L=256; o, 
L=384; A and L=512; O. 



With the objective of understanding how a single mutation spreads through 
genetic material, or how small hardware failures propagate on a computer ar- 
chitecture, we have compared two square lattices with identical sets of rules 
and almost identical Boolean variables at each site. The only difference be- 
tween them consists in a small modification of few central Boolean variables 
(typically less than 1% of the number of sites) . So one can observe the time- 
development of the two lattices and check the difference between them after 
a long time for each parameter value p . The difference between lattices is 
measured as the number of 'spins' that differ between the two lattices; this 
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difference is also known as Hamming distance d(t) between the lattice config- 
urations {(Tj(t)} and {pi(t)}: 



d(t) 



1 

N 



(2) 



Following earlier works [8,12,13], we call this difference the 'damage'. Further 
we define here the order parameter of the system as: 



ib = lim d(oo), 

d(0)— 

and one can also define a susceptibility as: 



X 



dd(oo) 
dd{0) 



(3) 



(4) 



It has been observed that ib as well as 1/x go to zero at some critical concen- 
tration p c in systems of dimensions greater than one for the short-range case 
Kauffman model similar to the para-ferromagnetic phase transition. In other 
words, for all p ^ p c , a small initial damage vanishes or remains small, i.e. be- 
longs to a small cluster of 'damaged spins', after a sufficiently long time. One 
says that the system is in the frozen phase. On the other hand, for all p > p c , 
a small initial damage spreads througout a considerable part of the system. 
Then one says that the system is in the chaotic phase. Of particular interest 
is, however, the border case p = p c where fractal properties appear (see next 
section). To know more about phase transitions in Kauffman's automata the 
reader should consult, for instance, Ref. [12]. As a way to find p c in both 
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Fig. 4. Log- log plot of ibL vs (p — p c )L x l v for the short-range case (q = 0.00). 
Lattice size L=256; o, L=384;A and L=512; O 

the short-range and the small- world case, we plotted the order parameter ib 
versus the parameter p for the case q = (short-range) and q = 0.05 and 0.10 
(small-world) for three different lattice sizes (L=256, 384 and 512) by flipping 
a small amount of regularly spaced central spins and taking an average over 
100 runs for each point. Everything was calculated for up to 3 x 10 3 time-steps 
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Fig. 5. Log- log plot of ipL vs 
(p - p c )L'^ u for the case (q = 0.05). 
Lattice size L=256; o, L=384;A and 
L=512; O 



Fig. 6. Log-log plot of tpL vs 
(p - PcJL -1 /" for the case (g = 0.10). 
Lattice size L=256; o, L=384;A and 
L=512; O 



on an Intel Pentium IV processor. The updating speed using a multi-spin cod- 
ing algorithm on an Intel machine which stores 64 spins in one computer word 
was about 2.44 updates per nanosecond for two configurations following the 
same iteration rules. Our results are shown in Figs. 2 and 3, repectively for 
q = and q = 0.05. The transition between the frozen phase and the chaotic 
phase is seen to occur around p c = 0.30 for the case q = and around p c = 0.17 
and p c = 0.16 for the cases q = 0.05 and q = 0.10, respectively. 

Table 1 





q = 0.00 


q = 0.05 


q = 0.10 


Pc 


0.30 


0.17 


0.16 


a 


0.31 ± 0.04 


0.78 ± 0.04 


0.74 ± 0.04 


7 


2.50 ±0.15 


4.75 ± 0.20 


4.65 ±0.15 



We estimate the critical exponents for each case by regarding the fraction 
of damaged sites (ip) on each side of the critical point after a sufficiently long 
time. In the chaotic range (p > p c ), this fraction goes to zero as (p — p c Y ■ 
While in the frozen range, the ratio of the final damage to the small initial 
damage (x) varies as (p c — p)" 1 . We have made a collapse of all data in order 
to check the consistency of the obtained results by making use of the following 
scaling law: 

iP(L,p) = L~^F((p-p c )L-^)i (5) 
where L is the linear dimension of the lattice and v is the exponent describing 
the divergence of the correlation length at p c . Fig. 4 shows the data collapse 
for the case q = at p c = 0.30 which is in good agreement with Ref. [14], 
whereas the Figs. 5 and 6 show the data collapse for the case q = 0.05 and 
q — 0.10 at p c = 0.17 and p c = 0.16 respectively. Indeed, from Figs. 4, 5 and 6 
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we observe that the points following the scaling relation given by Eq. 5 have 
a slope of straight line equal to critical exponent /3 for each case. This could 
be an indication of the quality of our data. We have evaluated 7 by means of 
the Eq. 4 at the transition threshold p c . Our estimates together with the p c 
values are shown in Table I. 



3 Fractal Dimensions 



In analogy with percolation theory, we have determined the fractal dimen- 
sions d act and d t both for the short-range and the small-world cases by calcu- 
lating the damage mass as well as the touching time on the lattice bounderies 
for several lattice sizes L at the critical point p c . 

The fractal dimension was defined by Stauffer [8] through the asymptotic pro- 
portionalities (L — > 00) in the following way: 

< M >~ L dact (6) 



<T t >~L d \ (7) 

where Eq. 6 says how the average damage increases with the lattice size L at 
the critical point, whereas Eq. 7 says how the average time for the damage 
spread increases with lattice size L at the critical point. The averages in Eqs. 
6 and 7 are taken over those runs that eventually touch the upper or lower 
boundary of the square lattices. In this paper as in Ref. [15] we flipped a 
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Fig. 7. Log-log plot of touching time 
vs L; •, damage vs L; o for the 
short-range case (q = 0.00) 




Fig. 8. Log-log plot of touching time 
vs L; •, damage vs L; o for the 
small-world case (q = 0.05) 



whole central line of spins at each time-step and let the system evolve until the 
final damage touches the lattice boundaries. Our results for the touching time 
T t and the damage mass M at the touching time ( averaged over 100 runs for 
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q = and over 1000 runs for q ^ 0) are shown as a log-log plot of time vs L 
in Figs. 7 and 8, respectively for the short-range (q = 0) and the small- world 
(q = 0.05) case. Our numerical results for the short-range case (q = 0) turned 
out to be d act = 1.79 and d t = 1.35. These results are in good agreement with 
Ref. [14] and obey the scaling relations: From table 1, we have (3 = 0.31 and 
7 = 2.48. Assuming hyperscaling one has dv = 7 + 2/3, and we get v = 1.55. 
Standard percolation theory, suggests d act = d — j3/v, predicting d act = 1.80. 
Our results for the small- world case (q = 0.05 and q = 0.10) were d act = 1.04 
and d t ~ 0. Therefore one can see easily that at q 7^ the system no longer 
obeys the scaling relations of percolation. 



4 Conclusion 

With the use of our algorithm based on the multi-spin coding technique, we 
simulate Kauffman's automata, analyzing how small central damages spread 
on square lattices both in the short-range case and the small-world case. In this 
work, we have seen that for the small- world case the speed of damage spreading 
is larger than for short-range case. We numerically evaluate critical exponents 
(3 and 7 for both cases. We also estimate the fractal dimension as being d act = 
1.79 for the short-range case and d act = 1.04 for small-world case. Where 
for the latter case we can no longer observe scaling relations for percolation. 
Finally, we expect that these and other more elaborated calculations will be 
helpful to understand more general problems concerning the propagation of 
simple defects in complex systems. 
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